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Abstract: We investigate the accretion process for different spherically symmetric space-time ge¬ 
ometries for a static fluid. We analyse this procedure using the most general black hole metric 
ansatz. After that, we examine the accretion process for specific spherically symmetric metrics 
obtaining the velocity of the sound during the process and the critical speed of the flow of the fluid 
around the black hole. In addition, we study the behaviour of the rate of change of the mass for 
each chosen metric for a barotropic fluid. 
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I. INTRODUCTION 

About two decades ara, astronomical observations 
such as CMB radiation [r Supernova Type la 0 and 
large scale structure data^, 01 showed that the Universe 
is undergoing an accelerating expansion period. This dis¬ 
covery was a revolution in cosmology and the agent re¬ 
sponsible for this effect was named dark energy. This 
energy has the strange property that produces repulsive 
gravitational effects and it violates the null and weak en¬ 
ergy condition 0, 0] . Surprisingly, observations suggests 
that around two-thirds of the total energy of the Uni¬ 
verse belongs to the dark energy Hi- However, the na¬ 
ture of this energy is not well-understood and nowadays 
is one of the most challenging problems in theoretical 
physics. Over the last two decades, a number of theorists 
have been trying to tackle this important problem. Sev¬ 
eral ideas have been proposed like cosmological constant, 
phantom energy, quintessence, k-essence, dynamic scalar 
fields, and others. Usually, dark energy is modelled us¬ 
ing a perfect fluid such as the pressure and energy density 
are related by a perfect fluid with a barotropic equation 
p = wp where w is the state parameter with w = —1 a 
cosmological constant, —1 < w < —1/3 for quintessence 
and w < — 1 for phantom models [^. 

Accretion is the process by which a massive astrophys¬ 
ical object such as a black hole or a star can take particles 
from a fluid from its vicinity which leads to increase in 
mass (and possibly angular momentum) of the accreting 
body 0. It is one of the most ubiquitous processes in 
the Universe. Indeed the stars and planets form only 
as a result of accretion in some inhomogeneous regions 
of gas and dust. Existence of supermassive black holes 
at the centres of giant elliptical and spiral galaxies sug¬ 
gests that such black holes could have evolved via accre¬ 
tion process. Other processes of formation of giant black 
holes such as merger of several small mass black holes 
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(or compact objects) or stellar collapse of several stars 
in a small domain leading to merger seem very remote 
possibilities. The most likely model of formation of mas¬ 
sive giant black holes is the accretion of dust or matter 
from nearby regions for sufficiently longer times. The for¬ 
mation of giant elongated astrophysical jets from active 
galaxies or small compact objects indicates the existence 
of large amounts of hot dust (likely to be in the form 
of accretion disk) around the regions of jet formation. 
However accretion process not always increases the mass 
of the compact source, sometimes the infalling matter is 
thrown away in the form of jets or cosmic rays. It is more 
probable that the accretion process may not be static and 
the velocity of free-fall and the energy density of the fluid 
change with time and position. The accretion of normal 
matter and dust onto compact objects is a well-studied 
problem, however the accretion of more exotic types of 
energy-matter is not so commonly probed including dark 
energy and stiff fluids. Since the Universe is dominated 
by dark energy, it is more pertinent to study the accretion 
of various forms of dark energy onto black holes. 


Bondi investigated accretion for compact objects using 
Newtonian gravity 0] . After the arrival of general relativ¬ 
ity, it became possible to investigate the accretion onto 
more compact objects such as neutron stars and black 
holes. Michel was the first who studied the relativistic 
accretion process for the Schwarzschild black hole M- 
Babichev et. al derived the fate of a stationary uncharged 
black hole in the phantom energy dominated universe 


Ilj . They concluded that the phantom energy falling 


onto a black hole will decrease its mass. We use the same 
formalism in order to study accretion process for different 
spherically symmetric spacetimes. More recently, Deb- 
nath presented a framework of static accretion onto 
general static spherically symmetric black holes, thereby 
generalising the work of Babichev et al [HI- We extend 
these previous studies using a more general ansatz for a 
static spherically symmetric spacetime which covers all 
possible static black hole (and non-black hole) solutions 
such as Schwarzschild-monopole, string cloud, Jannis- 
Newman-Winicour and dilaton (stringy charged) black 
hole. In particular, we investigate the parameters of fluid 
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and its effect on the black hole mass. 

This paper uses the geometrized units c = G = h = 1 
and the signature (—,+,+,+) for the metrics. This pa¬ 
per is organised as follows: In the second section, we de¬ 
rive a general formalism for the relativistic spherical ac¬ 
cretion process. In the third section, we study some spe¬ 
cific metrics such as black holes with topological defects, 
string cloud, JNW metric and a static charged black hole. 
For each case, we derive the velocity profiles, the speed 
of sound, the critical speed of the flow and the rate of 
change of the black hole mass for a barotropic fluid. We 
explain these parameters with the help of figures. Finally, 
we discuss the results of our paper. 

II. GENERAL FORMALISM 


lead to growth of black holes while the later ones lead 
to jets. Using the energy-momentum conservation law: 
0 = we find 

{p + p)u^^u^ + B{r)C{r) = A,, (5) 

Bir) 

where Ai is a an integration constant. 

By projecting the conservation law onto the 4-velocity 
= 0 , thereby contracting all indices, we can find 
the relativistic energy flux (or continuity) equation 

u^P,tJ. + {P + P)u^ ■,tj. = 0 . ( 6 ) 

In our case, we will assume that the pressure and energy 
density are related by a certain equation of state p = p{p). 
The last equation after simplification yields 


We consider the following metric ansatz for the general 
static spherically symmetric spacetime 


p' u' A' B' C 

p + p ^ ^ 


( 7 ) 


= —A{r)dt^ -\ - \—dr'^ + C{r){d0‘^ + sin^ 

B{r) 

where A{r) > 0, B{r) > 0 and C'(r) >0 are functions 
of r only. This metric ansatz with C(r) = was used 
by Debnath [ 13 . We are going to perform a general 
accretion formalism for those space-times geometries. 

Assuming the energy matter as a perfect fluid which 
is isotropic though inhomogeneous, the later assumption 
leads to variations in energy density at different positions 
near the accreting source. The energy-momentum tensor 
is given by 

Tfj.1, = ip + p)u^u^ +pgfj.,„ ( 2 ) 


Here, prime denotes differentiations with respect to r. 
Integration of (O yields 


uCir)^ p+p'(p) = -Aq, ( 8 ) 

where Aq is an integration constant, while a negative sign 
is introduced on the right hand side since u < 0 on the 
left hand side. Now, if we combine the above equation 
with (O we obtain 


ip + p)Vu^ + B{r)J 7 p+pip) = _^ = Ag, (9) 


where p is the energy density, p the pressure and is 
the four-velocity that in general will be 

u^ = ^ = {u\u\0,0), (3) 

where r is the proper time. We have that it® = 0 and 
= 0 due to the restriction of spherical symmetry. Note 
that all components of 4-velocity, pressure and energy 
density are functions of r only. Since the 4-velocity must 
satisfy the normalization condition it^it^ = — 1 we find 


( dt I + B 

dr V AB ’ 


( 4 ) 


where for simplicity we have named u = u'" = Due 
to the presence of square root, there are two possibil¬ 
ities: ii‘ > 0 (< 0 ), which respectively imply forward 
(backward) in time conditions. The former condition 
is necessary to preserve causality in the process: fluid 
must enter the black hole but not leave it. Moreover, 
to study accretion, we require it < 0 while for any out¬ 
ward flows It > 0. In the astrophysical context both in¬ 
ward and outward flows are important. The former ones 


where A 3 is a constant which depends on Ai and Aq. 
Due to spherical symmetry, we take 9 = tt/2 (fluid flow 
in the equatorial plane) which yields, y/—g = Cy/A/B. 
Furthermore, the equation of mass flux, 0 = = 

leads to 


where A 2 is an integration constant. If we divide Eq. ® 
with (fTOl) . we get another useful relation 


{p + p) 

P 



yju^ + B{r) 



( 11 ) 


where A 4 is another arbitrary constant. 

Taking differentials of Eqs. (HQl) and m and solving 
together, we obtain 



bJ u V A B J 



B' 

2{v? + B) 


dr = 0 , 


( 12 ) 
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where we have introduced the variable 


2 ^ d\n{p + p) 
d\np 


(13) 


It is important to mention that in Eq. (9) of Ref. [13 
a similar expression with C{r) = was derived which 
contained some erroneous signs. 

By taking the two brackets in (fT^ equal to zero, we can 
find the critical point of accretion located at r = r^. 
Thus, at the critical point we have 


C 


(K" - 1 ) 


- 


B{rc) - 
C'{r,) y2 
C{rc) ^ 


B{rc 


(14) 


B'{t,) 

2{ul + B{r,)) ■ 

(15) 


Here every function is evaluated at r = rc and Uc is 
the critical speed of the flow (velocity of the flow at the 
critical point). We can decouple the quantities and 
and obtain 




B{r)C{r)A' {r) 
2A{r)C'(r) ’ 
C{r)A’{r) 

C{r)A'{r) +2A(r)C'{r)' 


(16) 

(17) 


The speed of sound is found to be 


2 _ dp 

dp 


= Aa 


B{rc) 


A{rc)iul + B{rc)) 


-1. (18) 


Note that and cannot be negative, hence 


A'jrc) 

C'irc) 


> 0 . 


(19) 


From Eq. Cil), it follows that for given metric compo¬ 
nents, one can find the critical radius in the given space- 
time. 

The rate of change of the mass of the black hole M = 
dyrr^Tg, in the accretion process for a fluid around the 
black hole is given by [l^ 


Macc = 4^H3M2 (p + p), (20) 


where dot represents derivative with respect to time. 
This equation can be found by integrating the flux of 
the fluid over the surface of the black hole. We can see 
that the mass of the black hole will increase for any fluid 
such as p -I- p > 0 which accretes outside the black hole. 
On the other hand, if the fluid is a phantom dark en¬ 
ergy p + p < 0, then the mass of the black hole will 
decrease. It should be emphasized that this mass of the 
BH that appears in the corresponding metrics (see sub¬ 
sequent sections), where it is generally assumed that the 


mass is a pure constant. However, in realistic situations, 
mass of BH cannot remain fixed: accretion leads to in¬ 
crease in mass while Hawking radiation leads to the de¬ 
crease in mass. To connect the BH solutions with their 
astrophysics, one has to take into account the time depen¬ 
dence of BH mass. It is further assumed that this time 
dependence does not change the global geometry and the 
symmetries of the spacetime, therefore spacetime metrics 
remain static and spherically symmetric. If one is inter¬ 
ested to consider the combined effects of accretion and 
evaporation from BH simultaneously, then an additional 
term corresponding to Hawking radiation is added to the 
right hand side of Eq. (EOl), however it is not our concern 
at the moment. 


III. SPHERICALLY SYMMETRIC METRICS 
WITH HORIZONS 


In this section and subsequent, we are going to discuss 
the case in which Air) = B{r). All the spherically sym¬ 
metric metrics with horizons are included in this case. 
From Eqs. m and (HZl) we have 


u 


2 

c 


1^2 

c 


C{r)A'{r) 

2C'ir) ’ 

C{r)A'{r) 

C(r)A'{r) + 2A(r)C'{r) 


( 21 ) 

( 22 ) 


Although we are focused on the BH metrics with event 
horizons, the present analysis is not forbidden for hori¬ 
zonless spacetimes. In many cases, we are concerned with 
the properties of fluid flow (including critical velocity, 
sound speed in the fluid; evolution of fluid’s energy den¬ 
sity; speed of radial flow, etc.) and the change in mass of 
the accreting object and therefore the horizon does not 
directly involve anywhere. With this perspective, one can 
investigate the accretion phenomenon for more exotic ob¬ 
jects such as naked singularity and wormholes. Since the 
mass of BH is connected with the size of the horizon, the 
horizon expands due to accretion as well. 


A. Schwarzschild de Sitter Black Hole with a 
Topological Defect (Monopole) 

Dadhich et al. derived the metric of a Schwarzschild 
black hole with global monopole charge by relaxing 
asymptotic flatness and investigated several astrophysi- 
cal aspects [l3 • However in the presence of cosmological 
constant, the metric coefficients of this geometry become 
modify as follows [l3] 

A(r)=H(r) = 1- —- V, (23) 

r 3 

C{r) = (24) 

where rj and A are the parameter related to the scale 
of symmetry breaking and the cosmological constant re¬ 
spectively. This metric has a topological defect due 
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to the deficit solid angle around the scale of a gauge- 
symmetry breaking. Amani and Farahani [H investi¬ 
gated the phantom energy accretion onto Schwarzschild 
AdS black hole with global monopole and calculated the 
critical velocities of the flow. For this metric we find that 




3M - \rl 

6 rc 

\rl - 3M 
9M + 3Ar3 — 6 rc 


(25) 

(26) 


The condition (fT^ yields 


Arg - 3M 
3 {rf - 1 ) rg 

Using m, the speed of sound will be 


(27) 


cl = Aa 


2 rr 


2rr — 3M — Arg 


- 1 . 


(28) 


It is possible to integrate the conservation laws and ob¬ 
tain analytical expressions of the physical parameters. 
Here, for simplicity we will again study the barotropic 
case where the fluid has an equation of state such as 
p(r) = wp{r). Using ([5]), Q and (fTTl) we can find di¬ 
rectly one set of solutions 


p(r) 

p(r) 


u{r) 


\/%AqA 2 

{rj^ — 1) r 3 / 2 -\/Ag ( 6 M -I- Ar^ — 3 r) -|- 3 A|A|r 
%/3Ao(Ao — A 2 ) 

{rf — 1) ( 6 M + Ar® — 3 r) -|- 3 A|A|r 

y/Ag ( 6 M + Ar 3 - 3 r) -t 3 A^A|r 
\/3Ao-v/r 


(29) 

(30) 

(31) 


The other set of solutions are unphysical solution be¬ 
cause they have p < 0 for w > — 1. Now, if we only use 
and CD and we assume a barotropic fluid such as 
p = wp we obtain 


p(r) = 


V3A2 (w + 1 ) 


( 7 ]^ — 1) sjr ( 3 A 4 3)) GM{w + 


(32) 

x/r (3AJ -I- (UJ -I- 1)2 (Ar2 - 3)) -I- 6M{w + 1)2 
= -■ (33) 

V3x/r{w + 1 ) 


Note that if we compare both solutions, we can notice 
that if we replace Aq = 1 -I- w and A 2 = 1 in (IMli-dsD we 
obtain the same solution as the equations above. There¬ 
fore, if we choose Aq = l-|-w and A 2 = 1 there is an equiv¬ 
alence between solving the system of equations given by 
® , (HI and (fTTl) and solving the system of equation given 
by ®, (CD and p = wp. Thus, in order to study the ac¬ 
cretion process for different kind of barotropic fluids, we 
will choose Aq = 1 -I- w and A 2 = 1. Using this argu¬ 
ment, for the definitions of the constants we must have 
that A 3 = A 2 A 4/(1 -I- ic). This interesting fact is be¬ 
cause (O is just a contraction of ([51), and therefore there 
are not more degrees of freedom in our equations and it 


seems like the only possibly solution for the system is 
to have a barotropic fluid (with different possibles state 
parameters). For all the following subsections, it can be 
showed that there is also an equivalence too. 

Fig. (P) represents the absolute value of the velocity pro¬ 
file for different values of the state parameter. Here 
w < —1, —1 < u> < —1/3, w = —1, w = 0, w = 1 refer 
to phantom energy, quintessence, cosmological constant, 
dust and stiff matter, respectively. The dots denote the 
critical velocity of the fluid at a certain critical radius. 
The behaviour shown in these figures is consistent with 
previous results reported in [I^, 171. It can be observed 
that when w > 0 increases, the location of critical radius 
shifts to the left. Thus the infalling fluid acquires super¬ 
sonic speeds closer to the black hole. We also took A < 0 
in figures, which correspond to anti-de Sitter background. 
Fig (P represents the behaviour of the energy density of 
fluid for different values of the state parameter. It can 
be seen that the energy density of the fluid decreases as 
it approaches the black hole when w = —2 and w = —1.5 
while it increases in the remaining cases. In the former 
cases, the energy density becomes negative which violate 
the energy conditions too. However, in absolute value, 
the energy density always increases for all cases. Asymp¬ 
totically the fluid energy density approaches zero at infin¬ 
ity while it approaches to a maximum near the black hole 
due to the strong (and perhaps quantum) gravitational 
interaction. 



FIG. 1: Velocity profile versus x = -^ for A 4 = 10, A = —0.73, 
M = 1 and different values of the state parameter. For each 
figure, the critical radius are marked with a point for every 
chosen state parameter. 


From (I20I) . the rate of change of the mass of the black 
hole due to the accretion process for a barotropic fluid 
becomes 


M 


_ 4x/3'KAyA4M^(w + l) _ 

(r;2-l)r3/2y/r(3A|-|-(iu-|-l)2(Ar2_3))-|-6M(u;-|-l)2 


.(34) 
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FIG. 2: Energy density versus x = -^ for A 2 = 1, A 4 = 10, 
A = —0.73, M = 1, r] = 0.5 and for different values of the 
state parameter. 


Fig (l3|) represents the rate of the change of mass of 
the black hole for different values of the state parameter 
and same initial black hole mass. For small values of x, 
near the black hole, M > 0 showing that mass of black 
hole increases for w = 0 and w = —0.5 being matter 
and quintessence respectively. In other words, mass of 
the black hole will increase for matter and quintessence 
accretions. However for stiff matter accretion, M > 0 for 
both small and large values of x, thereby increasing in 
mass of black hole. It should be noted that quintessence, 
dust and stiff matter satisfy the relativistic energy condi¬ 
tions for energy-matter which ensure the increase in BH 
mass and satisfy the second law of black hole thermody¬ 
namics. It is also apparent that M < 0 for phantom-like 
equations of state such as w = —1.5,—2, thus mass of 
BH will decrease by the accretion of phantom like fluids. 



FIG. 3: Rate of change of the mass of the black hole versus 
X = -^ for a barotropic fluid for different values of the state 
parameter, A 2 = 1, H 4 = 10, M = 1, A = —0.73 and rj = 0.5. 


It is important to remark that we have another possi¬ 
bility to study the accretion process for this metric, which 
is by taking H 2 < 0 and using (EH) and (1331) with a minus 
of difference. In that case we have that p > 0 for w > — 1 
as we required but u{r) < 0 for w < — 1. However, the 


case H 2 < 0 seems to be an unnatural choice since we 
can notice that we will need a barotropic case such as 
p = —wp to have identical solutions to solving ([H, dH 
and CH) directly and setting Aq = 1 + w and A 2 = — 1. 
In other words, it is the same to have Aq = 1 + w and 
H 2 = — 1 and solving the system of equations (O, dH and 
m than solving the system of equations dH-CT and 
p = —wp. Therefore, we will have the same behaviour as 
we discussed above since for y 42 = — 1, the state param¬ 
eter w will be minus the state parameter in the case of 
A 2 = 1. Thus, for all the following subsections, we are 
going to find the solutions using ®,dlll) and p = wp and 
we will take A 2 = 1 . 


B. Schwarzschild Black Hole in a String Cloud 
Background 


The metric for a string cloud arises from the Nambu- 
Goto action of a string evolving in spacetime. In 4- 
dimensions, the metric of a Schwarzchild black hole in 
a string cloud background is given by [13 

A{r) = B{r) = 1 — — a, (35) 

r 

Or) = (36) 


where a is the string cloud parameter which is a 7 ^ I and 
M is the mass of black hole and it appears as a constant 
of integration during the solution of field equations [l 8 l| . 
The limit a = 0 gives the Schwarzschild solution. The 
radial velocity and the energy density for a barotropic 
fluid are given by 


u{r) = 

p{r) = 


\/r {Al -I- (a - l)(ic -b 1)^) -b 2M(w + 1)^ 


\/r{w + 1 ) 

A2{w + 1 ) 


(37) 


{A\ -I- (a — l)(w + 1)^) -b 2M{w 1)^ 

(38) 


Exactly as before, we have two set of solutions for u{r) 
and p{r). For this case, the critical values are given by 

3M(w -b 1)2 

“ ~2{Al + {a-l){w + lYy 
, M 


° 3M-b2(a-l)rc' 

As a special case, if A 4 = 1-l-w, then the above equations 
yield Tc = —3M/l-|-2a, which for a = 0 gives the critical 
radius for the Schwarzschild black hole. Since is a 
function of w, its value will be different for every fluid. 


(39) 

(40) 

(41) 
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FIG. 4: velocity profile versus x = -^ for A 4 , = 1, a = —5, 
M = 1 and different values of the state parameter. For each 
figure, the critical radius are marked with a point for every 
chosen state parameter. 



FIG. 6 : Rate of change of the mass of the black hole versus 
X = for a barotropic fluid for different values of the 
state parameter, A 2 = 1, A 4 = 1, M = 1 and a = —5. 
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FIG. 5: Energy density versus a; = for ^42 = 1, A 4 — 1, 
a = —5, M = 1 and for different values of the state parameter. 


From Fig. o, it can be seen that the radial velocity of 
the fluid is negative for the cases w = —1.5 and w = —2 
while it is positive for w = —0.5, 0,1. If the fluid flow is 
inwards then u > 0 is not allowed and vice versa. In the 
case of quintessence, the fluid velocity is zero (at rest) at 
X = 1 while it attains critical velocity at approximately 
X = 0.7. Similarly for dust and stiff matter, the critical 
values of velocity occur at approximately x = 0.3. 
Similarly, accretion of phantom fluids lead to critical 
flows after x = 0.7. In Fig. dSD, we plot the evolution of 
energy density of fluids in the vicinity of the black hole. 
It is apparent that for fluids satisfying the weak and 
dominant energy conditions such as dust, quintessence 
and stiff fluid, the energy density increases as the fluid 
moves towards the BH while the reverse happens for 
w = —1.5 and w = —2. In Fig. we make a plot of 
rate of change of mass of BH against x and find that 
mass of the BH will increase due to the fluids satisfying 
energy conditions while it decreases for fluids violating 
the same energy conditions. 



Eq. m is completely satisfied since M > 0. The 
speed of sound is given by 


C. Janis-Newman-Winicour Spacetime 


= Aa 



- 1 . 


(42) 


Thus, the rate of change of the mass of the fluid is 


M = 


A-kAIA4,M'^{w -I- 1) 


(Af + (a — l)(w -|- 1)^) + 2M{w + 1)^ 


(43) 


The physical existence of naked singularities is still 
questionable, however there are few exact solutions of 
Einstein field equations that represent naked singular- 
ity. These include the JNW, JMN and gamma metrics 
[1^ . The Janis-Newman-Winicour (JNW) solution is ob¬ 
tained as an extension of the Schwarzschild space-time 
when a massless scalar field (with vanishing potential) is 
introduced. Thus this solution is not a vacuum solution 
. Here the coordinate singularity in the Schwarzschild 
space-time becomes a naked singularity in JNW space- 
time. Joshi and collaborators (see also references 
therein) have investigated several astrophysical features 
distinguishing a black hole from a naked singularity. The 
gravitational leasing features due to naked singularity 
have also been explored [ 2 l|. 
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There are various forms of JNW metric available in 
literature, however we adopt the following form 


A{r) = B{r) = 

-2r - ro[fi - l)l i 

-2r-1-ro(Ai + 1)J 

(44) 

C{r) = 

1 [2r-f ro(Ai + l)]^+^ 

(45) 

4 [2r - ro(Ai - 1)]^* ^ 


where /i is the scalar charge which is defined in the inter¬ 
val (1, oo) and ro is related to the mass by ro = 2M. The 
Schwartzschild metric can be obtained by taking = 1 
after a coordinate transformation. The scalar field is 


2 ^ - roifJ. - 1) 

H 2r + ro{fi + l) 


where a and fi are related by 


— 1 -f 


327ra 


2 


(46) 


(47) 


In this case, the 4-velocity and the energy density of the 
barotropic fluid are 


- , \/A\ — {w + l)2(2r — pro -I- ro)i/'"(2r + pro -|- 

=-7m- 


(48) 



0.50 0.55 0.60 0.65 

X 


FIG. 7: Velocity profile versus cc = ^ for A 4 = 0.3, 
M = I,^ = 2 , ro = l and different values of the state 
parameter. For each figure, the critical radius are marked 
with a point for every chosen state parameter. 


critical velocities for different values of x, however in¬ 
terestingly we do observe some symmetry in the profile 
of the velocity curves and the location of critical points. 
After passing the critical point, the fluid flow becomes 
super-sonic or trans-sonic. In Fig. ([ 8 ]), the graph of en¬ 
ergy density of different fluids in the neighbourhood of 
naked singularity. The fluid energy density for dust and 
quintessence rises in the near vicinity of singularity and 
goes to zero asymptotically. Similarly, we also see that 
energy density of dust and quintessence remains positive 
while for phantom fluids, it is negative. 


P(x) 



FIG. 8 : Energy density versus x = lor A 2 = 1, A 4 = 0.3, 
p = 2 , ro = 1 and for different values of the state parameter. 


In addition, for this metric we find that 

(50) 

(51) 

1.(52) 

Using the condition (fT^ we obtain 


U" = 


tq f 2(2re + rp) 

4rc V2rc + (a* + l)?’o 
rp 

4rc -I- ro ’ 

2 A 4 . 


- 1 


1//X 


2rc(l - /r)ro \ 


4rc -I- ro V2rc + (a* + l)?'o 


2ro(2rc - A^ro + ro)e' 


i_i f 2rc - fJ-ro +ro\ 


1 /m 


2rc -I- fJ.ro + ro / 
(2rc -I- fjro + roY 


> 0 . 

(53) 


For this case, the fluid will have a change of its mass as 
follows 


In Fig. we plot the radial velocity of the fluid 
for different values of the state parameter. The fluids 
have zero radial velocity in the asymptotic limit, how¬ 
ever they acquire non-vanishing velocities as they ap¬ 
proach the naked singularity. Due to strong gravitational 
attraction near the naked singularity, the fluids achieve 


00 — — 1 — 

^ _ 167tA2A4M^{w + l)(2r + (1 - m)^o) ^ (2r + (/x + l)ro) ^ 

Y- (w + l)2(2r + (1 - /r)ro)i/f‘(2r + (m -I- 

(54) 

From Fig. ( 0 ), that the mass of the JNW singularity 
increases in the cases of dust and quintessence while it 
decreases for phantom fluids. 





























FIG. 9: Rate of change of the mass of the black hole versus 
X = ^ for a barotropic fluid for different values of the state 
parameter, A 2 = 1, A 4 = 0.3, M = 1, fi = 2 and ro = 1. 


D. Charged Black Hole in String Theory 


Gibbons and Maeda [ 2 ^ and later on independently 
Garfinkle, Horowitz, and Strominger [23] discovered a 
static spherical symmetric charged black hole in low en¬ 
ergy effective theory of heterotic string theory in four 
dimensions. Sharif and Abbas [l^l investigated the dy¬ 
namical behavior of phantom energy w < —1 near stringy 
charged black hole. They deduced that this mechanism 
reduces the mass of the black hole, though the black 
hole cannot be convertible to the extremal black hole. 
We again analyze this black hole for accretion dynamics 
without restricting w. We write down the stringy charged 
metric in the following form [^ 


B{r) 

2M 

r 

(55) 

Cir) 

_ A Q^exp(-2$o)\ 2 

V Mr / ’ 

(56) 


where Q is the charge parameter and $0 the asymptotic 
value of the dilatonic field which can be set to zero under 
certain cases [l^j. In this case, we have 


u 


2 

C 




C 


2 

s 


MQ"^ - 

Q'^r — 2Mr^e^'^° ’ 

M{Q'^ - Mre‘^^«) 

Q2(r -M) + Mre2-J>o(3M - 2r) ’ 


r (2Mre2‘fo - Q^) 


Q‘^{M -r) + Mre'^’^o (2r - 3M) 


(57) 

(58) 

1.(59) 


The radial component of the 4-velocity of the fluid and 
the energy density are respectively 


u(r) = — 


\/r (A| — (w + 1)2) -I- 2M{w + 1)2 


^(w + 1 ) 


(60) 


p(r) = - 


A2Me^‘^«iw + l) 


^\/r (A| — (u) -|- 1)2) -I- 2M{w + 1)2 [Mre^^o — 

(61) 


In Figs. (fTOl) and dill) , we plotted the radial component 
of the fluid flow velocity against the dimensionless radial 
coordinate x. We also represent the location of critical 
velocity by thick dots whose presence in the fluid flow 
profiles show that critical flow is a generic property of 
all fluids in the strong gravitational regimes. All velocity 
prohles u{x) > 0 are physically forbidden since nothing 
can escape from the black hole. However flows of dust, 
stiff matter and quintessence is permissible under our 
assumptions of inward flows m < 0 for all x. 



FIG. 10: Velocity profile versus x = ^ for A 4 = 2.1, M = 1, 
Q = 1.5, $0 = —0.1 and different values of the state param¬ 
eter. For each figure, the critical radius are marked with a 
point for every chosen state parameter. 



FIG. 11: Velocity profile versus x = -^ for A 4 = 2.1, M = 1, 
Q = 1.5, 4>o = —0.1 and for w = 1.0. The critical radius are 
marked with a point. 


In Figs. (na) and ca), we plotted the energy den¬ 
sity profiles of the fluids in the vicinity black hole. For 
phantom-like fluids, the energy density is negative while 
for the quintessence, dust and stiff fluids the energy den¬ 
sity is positive. It also shows that fluids will have higher 
densities compared to their asymptotic values. 
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X 

FIG. 12: Energy density versus x = ^ ior At = 2.1, A 2 = 1, 
M = 1, Q = 1.5, 'Fo = —0.1 and for different values of the 
state parameter. 



FIG. 14: Rate of change of the mass of the black hole versus 
X = ^ for a barotropic fluid with different values of the state 
parameter and the following parameters: A 4 — 2.1, A 2 = 1, 
M = 1, Q — 1.5 and "Fo = —0.1. 


P(x) 




FIG. 13: Energy density versus x = -^ for A 4 — 2.1, A 2 — 1, 
M — 1, Q — 1.5, Fq = —0.1 and for ui = 1. 


By using (1191) we find the following constraint for the 
critical radius 


FIG. 15: Rate of change of the mass of the black hole versus 
X = -^ for a barotropic fluid with w = 1 and the following 
parameters: A 4 = 2.1, A 2 = 1, M = 1, Q = 1.5 and Fq = 
- 0 . 1 . 


IV. DISCUSSION 


2 ^/2g2<I.o 

r2 (2Mre2^o - Q^) ^ 

Therefore, the rate of change of the mass is given by 
A-KA'iA 4 M^e'^^^ (w + 1 ) 

M =- ^ ^ ’ -. 

\/r\/r(A4 — w — 1)(A4 + + 1) + 2M{w + 1 )^ (Mre^^o — 

(63) 

In Figs. (HI and (USD, the rate of change in BH mass M 
is plotted against x for different values of state parameter. 
In all cases, it is observed that the BH mass will increase 
by the accretion of quintessence, dust or the stiff fluid 
when X is small. Strikingly, the maximum rate of increase 
in BH mass arises due to quintessence followed by dust 
and stiff fluid. 


We have proposed a most general framework for the 
study of spherical accretion onto spherically symmetric 
compact objects. We explored various features of the 
fluid flow near black holes. Although we assumed all the 
fluids to satisfy a certain linear equation of state with 
distinct state parameters, the above analysis can alter¬ 
natively be performed without specifying the equation 
of state, since the number of equations are sufficient to 
close the system of equations. Interestingly we found a 
link between the barotropic equation of state and the 
conservation laws. However the equation of state helps 
us identify which kind of fluid is falling onto the black 
hole. In other words there is no multiple accretion sce¬ 
nario into play, though it can be done. In particular, we 
focused on dust, stiff matter, quintessence and phantom 
dark energy accretion onto black holes unlike previous 
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studies in literature which focused on only one kind of 
test fluid. We did not consider the case of cosmological 
constant or vacuum energy since its accretion does not 
alter the evolution of black holes. It is found that dif¬ 
ferent fluids with distinct state parameters have different 
evolutions in the black hole backgrounds. Certain flu¬ 
ids acquire positive or negative energy density near the 
black hole while some fluids cause black hole masses to 
increase or decrease. Although we plotted all the fluid 
behaviours via single graphs, it is assumed that a sin¬ 


gle test fluid is accreted at one time. For a future work 
it will be interesting to perform similar analysis with a 
dynamically spherically symmetric geometry or with a 
non-static fluid. 
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